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ABSTRACT 

The  Office  of  Naval  Research  (Code  222)  is  sponsoring  an  experi- 
mental program  in  the  Northeast  Pacific  Ocean  to  investigate  long-term 
variations  in  propagation  loss  and  travel  time  between  a fixed  source 
and  fixed  receiver  over  extremely  long  paths.  This  report  considers 
environmental  factors  which  are  expected  to  be  important  for  acoustic 
modeling  and  prediction,  particular  attention  being  given  to  mesoscale 
features,  or  "ocean  weather."  Existing  and  planned  environmental  data 
sources  and  formats  are  reviewed  in  terms  of  the  requirements.  Cost- 
effectiveness  is  emphasized.  Monthly  weeklong  shallow  XBT  (480  m) 
sampling  of  the  3835-km  Oahu  to  San  Francisco  path  by  the  National 
Marine  Fisheries  Service  (NOAA)  ship-of-opportunity  program  is  found 
to  serve  well  as  a basis.  However,  the  sampling  depth  must  be  increased 
beyond  the  sound  channel  axis  at  about  800  m,  and  the  spacing  of  drops 
should  be  reduced  to  about  30  km  in  the  vicinity  of  sharp  temperature 
and  salinity  gradients.  A program  has  been  implemented  to  merge  the 
taped  historical  oceanographic  station  data  file  (NAVOCEANO)  with 
synoptic  XBT  and  salinity  data  (NMFS) . Available  bathymetric  (NORDA) 
and  Nh  tidal  (NSWC)  data  banks  appear  to  be  quite  satisfactory.  It  is 
highly  desirable  to  obtain  long-term  time-series  temperature  variability 
data  down  to  depths  below  the  sound  channel  axis  using  thermistor  arrays 
at  a suitable  location  along  the  path.  The  use  of  satellites  and  remote 
sensing  aircraft  to  supplement  the  existing  data  acquisition  programs  is 
strongly  recommended.  Existing  computer  models  of  extremely  long-range 
acoustic  transmission  loss  verify  relatively  large  measured  range- 
dependent  effects  upon  crossing  fronts,  currents,  and  eddies.  The 
results  of  previous  long-term  time-series  measurements  of  propagation 
loss  for  the  1250-km  Eleuthera-Bermuda  fixed  path  at  406  Hz  are  pre- 
sented as  a reference. 
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I.  INTRODUCTION 

The  Office  of  Naval  Research  (Code  222)  is  sponsoring  an  experimen- 
tal program  in  the  Northeast  Pacific  Ocean  to  investigate  long-term  vari- 
ations in  propagation  loss  and  travel  time  between  a fixed  source  and 
fixed  receiver  over  extremely  long  paths.  Particular  attention  is  given 
to  a fixed  acoustic  path  between  Oahu,  Hawaii,  and  San  Francisco  (3835  km). 
The  purpose  of  this  ONR  study  is  to  determine  the  relationship  of  these 
variations  to  large  space-time  scale  ocean  processes.  The  acoustic  data 
should  also  be  useful  for  the  study  of  fluctuations  and  variations  re- 
sulting from  other  ocean  effects,  e.g.,  internal  waves,  tides,  finestruc- 
ture,  surface  waves  and  currents,  current  gradients,  and  the  presence  of 
bottom  topography. 

This  report  considers  the  environmental  factors  that  are  expected 
to  be  significant  for  studying  long-term  acoustic  variability.  The  pur- 
pose of  the  present  work  is  to  recommend  a cost-effective  environmental 
data  acquisition  program  taking  into  account  the  most  significant  ocean 
variables  that  can  be  applied  to  acoustic  modeling  in  support  of  the  ONR 
fixed-source  fixed-receiver  propagation  measurements.  The  focus  of  the 
work  has  been  to  uncover  existing  data  sources  and  data  banks  and  to 
examine  the  adequacy  of  existing  measurement  and  analysis  capabilities 
which  could  be  applied  to  the  problem.  The  first  stage  has  been  planned 
and  is  being  implemented  by  ONR  in  collaboration  with  the  National  Marine 
Fisheries  Service  (NMFS,  NOAA) , the  Naval  Oceanographic  Office  (NAVOCEANO) , 
and  the  Fleet  Numerical  Weather  Central  (FNWC) . Work  is  continuing  on 
expanding  the  Northeast  Pacific  areal  coverage,  and  on  relating  mesoscale 
surface  features  to  features  at  depth. 

Statistical,  spectral,  dynamic  structural,  and  historical  data  meth- 
ods are  being  applied  to  the  problem  of  describing  and  forecasting  "ocean 
weather"  for  the  present  acoustical  purpose.  The  term  "ocean  weather"  is 
appropriate  because  the  time  and  length  scales  that  are  associated  with 
the  warm  and  cold  mesoscale  ocean  features  are  reminiscent  of  cyclones 
and  anticyclones  observed  in  atmospheric  weather  systems.  The  material 
on  water  masses,  fronts,  currents,  tides,  bottom  topography,  and  sound 
speed  provinces  summarizes  data  and  data  systems  presently  available, 
albeit  difficult  to  obtain.  Next  covered  are  programs  in  progress  in- 
volving sampling  of  the  internal  temperature,  density  and  sound  speed 
field,  dynamic  height,  currents,  and  wave  heights  related  to  instrumen- 
tation and  observational  platforms  such  as  ships,  buoys,  aircraft  and 
satellites.  Finally,  presently  available  experimental  acoustic  data 
over  long-range  paths  in  other  parts  of  the  world  oceans  are  introduced, 
along  with  performance  estimates  provided  by  existing  acoustic  models 
for  fixed  and  moving  sources,  narrow  and  broad  band. 
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II.  NORTH  PACIFIC  OCEAN  WATER  MASSES,  FRONTS  AND  CURRENTS 
I I -1 . Schematic  Map  of  Main  North  Pacific  Fronts 

First  let  us  review  the  oceanographic  features  of  the  North  Pacific 
Ocean  and  how  these  relate  to  the  problem  of  monitoring  the  environment 
of  an  acoustic  path  between  Oahu  and  San  Francisco.  We  are  particularly 
interested  in  the  time  and  length  scales  of  these  environmental  features. 

Figure  I I - 1 shows  the  basic  frontal,  water  mass,  and  water  current 
regions  of  the  North  Pacific.  The  region  between  the  Subtropical  Front 
and  the  Subarctic  Front  is  called  the  Transition  Zone.  The  acoustic  path 
starts  at  Oahu  in  subtropical  water,  known  as  Eastern  North  Pacific  Central 
Water,  crosses  the  Subtropical  Front  near  the  position  of  Ocean  Weather 
Ship  (OWS)  NOVEMBER  at  about  the  halfway  point  (30°N,  140°W)  , passes 
through  a section  of  the  Transition  Zone,  enters  the  California  Current 
region,  and  ends  in  the  coastal  zone  of  San  Francisco. 

Table  1 1 - 1 provides  a comparison  of  the  horizontal  acoustic  and 
oceanographic  properties  of  the  main  North  Pacific  frontal  and  current 
zones. 


Table  II-l.  Characteristics  of  North  Pacific  fronts,  based  on  horizon- 
tal sampling  at  30  km  intervals.  All  gradients  are  ex- 
pressed in  terms  of  the  Nyquist  sampling  interval.  (from 
Ref.  II-l;  reproduced  with  permission) 

North  Pacific  Fronts 


Characteristics* 

Kuroshio 

Oyashio 

Subarctic 

Subtropical 

Do Id rum 

Temperature  gradient  (“C  per  60  km) 

6 

9 

8 

4 

1 

Salinity  gradient  (i.  per  60  km) 

Sound  velocity  gradient  (m  s'*  per  b0  km) 

0.6 

1.5 

1.2 

0.5 

1.0 

24 

39 

28 

12 

1 

Density  gradient  (kg  m'^  per  60  km) 
Baroclinic  current  (m  s'*) 

0.8 

0.2 

0.2 

0.8 

0.7 

0.6 

0.2 

0.4 

0.5 

0.5 

Baroclinic  shear  (s'*) 

2 x 10'5 

8 x 10'6 

7 x 10'6 

10'5 

2 x 10': 

'Maximum  values 


Figure  11-2  is  taken  from  The  Oceans,  by  Sverdrup,  Johnson  and 
Fleming,  and  shows  the  distribution  of  water  masses  as  identified  by 
their  T-S  curves.  We  see  that  the  acoustic  path  goes  through  Eastern 
North  Pacific  Central  Water,  the  Transition  Region,  and  the  California 
Current . 
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MAIN  OCEANIC  FRONTS 


directions . Letters  refer  to  fronts:  A,  Kuroshio  front;  B,  Oyashio  front; 

C,  subarctic  front;  D,  subtropical  front;  E,  doldrum  front.  (from  Ref.  II-l; 
reproduced  with  permission) 


and  temperature  distributions , as  denoted  by  their  T-S  curves.  (from  Ref.  II-2; 
reproduced  with  permission) 
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11-2.  Fixed  Acoustic  Path  Traverses  Ocean  Sub-Gyre 

An  important  property  of  the  Oahu  to  San  Francisco  acoustic  path  can 
be  seen  in  Figure  II-3.  For  the  most  part,  the  path  passes  through  an 
oceanic  sub-gyre  formed  by  the  splitting  of  the  primary  clockwise  circula- 
tion pattern:  easterly  trade  winds,  Kuroshio,  westerlies,  California  cur- 
rent, and  into  two  wind-driven  gyres.  The  split  occurs  just  to  the  west  of 
the  Hawaiian  Islands.  The  boundary  of  the  gyre  migrates  somewhat  over  the 
year,  and  from  year  to  year,  depending  on  the  position  of  the  atmospheric 
semipermanent  high  pressure  region  that  overlies  the  Eastern  North  Pacific. 

NAVOCEANO  SP  1402-NP8,  "Surface  Currents--North  Central  North  Pacific 
Ocean"  (Ref.  1 1 -4 ) , and  others  in  the  same  series  provide  a summary  of 
historical  surface  current  vectors  by  month  for  1°  x 1°  squares.  Such 
data  may  be  useful  for  studying  the  effects  of  water  motion  and  circula- 
tion on  acoustic  propagation. 

Figure  11-4  shows  the  transection  of  three  radials  from  Oahu  to 
Los  Angeles,  San  Francisco  and  Seattle  through  the  same  water  mass  regions. 
The  radial  to  San  Francisco  seems  to  be  a turning  locus  for  the  cast-west 
section  of  the  Transition  Zone  to  a more  north-south  orientation.  The 
observations  were  taken  as  part  of  the  ship-of-opportunity  XBT  measure- 
ment program  of  the  National  Marine  Fisheries  Service  (NOAA) . 
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160°  140°  120°  W 


Figure  II-4.  Sections  of  expendable  bathythermograph  observations  were 
taken  on,  or  near,  the  great  circle  route  between  Oahu 
and  San  Francisco.  The  surface  water  masses  in  the  region 
are  shown  schematically . (from  Ref.  II-5) 
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III.  TIDES  OF  THE  NORTHEAST  PACIFIC 
I I I - 1 . Introduction 

Tidal  effects  are  very  important  for  the  fixed  acoustic  path:  They 
stretch  or  constrict  vertical  ocean  cross  sections,  cause  the  propagation 
of  internal  waves,  and  cause  currents  to  flow  in  response  to  their  forces 
Figure  1 1 1 - 1 shows  semidiurnal,  or  M2 , cotidal  lines  of  the  world  ocean 
as  drawn  in  1904.  This  is  the  most  important  tidal  component.  The 
cotidal  lines  are  lines  of  constant  tidal  phase.  The  points  from  which 
they  radiate  are  called  amphidromic  points  or  amphidromes.  As  the  earth 
turns,  the  cotidal  lines  appear  to  rotate  about  the  amphidromes. 

1 1 1 - 2 . Schwiderski  Model  for  M2  Tide 


The  most  complete  predictive  study  of  global  ocean  tides  has  been 
given  by  Dr.  E.W.  Schwiderski  (Ref.  III-l).  He  has  been  able  to  develop 
and  test  a discrete  tide  model  for  a 1°  x 1°  graded  grid  system  in 
connection  with  a hydrodynamically  defined  bathymetry.  The  Laplace 
tidal  equations  are  augmented  by  turbulent  friction  terms  with  novel 
mesh  area  (latitude  and  depth)  dependent  eddy-viscosity  and  bottom- 
friction  coefficients.  The  well-known  astronomical  tide-generating 
forces  are  modified  by  effects  due  to  solid  earth  tides  and  ocean  tidal 
loading. 

He  has  computed  the  M2  tide  in  map  form  for  four  typical  30°  (long.) 
by  50°  (lat.)  ocean  areas.  Tables  1 1 1 - 1 and  1 1 1 -2  show  the  amplitude 
and  phase  for  the  Northeast  Pacific.  It  is  estimated  that  the  tabulated 
tidal  charts  permit  a prediction  of  the  M2-tide  elevation  of  the  ocean 
surface  over  the  geoidal  level  with  an  accuracy  of  better  than  5 cm 
anywhere  in  the  open  ocean  and  with  somewhat  less  accuracy  near  rough 
shorelines.  With  the  forthcoming  construction  of  the  lesser  tidal 
constituents,  the  total  tidal  prediction  error  can  be  kept  below  a 10  cm 
bound.  The  lesser  tidal  components  to  be  modeled  next  are  S2,  , and  0^ . 

1 1 1 - 3 . Tidal  Data  Checks  of  Schwiderski  Model  in  the  Northeast  Pacific 


Two  independent  checks  of  Schwiderski ' s model  against  known  measure- 
ments in  the  Northeast  Pacific  show  remarkable  agreement. 

(a)  Larsen  and  Irish  (Ref.  II I -la)  made  measurements  of  the 
M2  tide  at  Cobb  Seamount  (46°  46.4'N;  130°  48.8'W).  The 
comparison  of  their  data  with  Schwiderski ' s prediction  is 
shown  in  Table  1 1 1 -3 . 
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KEY  TO  TABLES  111-1  AND  I II -2 


M * Longitude  East  (*) 

N * Colatitude  (d) 

£ * Amplitude  (cm)  * half  range  between  high  and  low  water. 

5 * Phase  (°)  « arrival  time  (1°  * 2.070  min.)  of  tidal  crest 
after  the  moon's  passage  over  Greenwich  meridian. 

$ * Circled  cross  marks  amphidromic  point  of  zero  amplitude. 

Subbars  mark  empirical  tide  input  values. 

i i Subbrackets  mark  observed  near-shore  deep-sea  input  data. 

. . . . Wavy  underlines  mark  off-shore  deep-sea  tide  gauge  stations 
with  observations  listed  in  Tables  4a  and  b (see  also  Tables 
3 and  2).  [tables  in  Ref.  III-l) 

Estimated  accuracy:  about  S cm  in  open  ocean. 


Table  III-3.  tide  component  at  Cobb  Seamount. 


Amplitude 

Phase 


Larsen  and  Irish 
(measured) 

81 .07  cm 

242.0° 


Schwiderski 

(predicted) 


82.1  cm 


245.0° 


(b)  Figure  1 1 1-2  shows  the  M2  cotidal  chart  from  Munk, 
Snodgrass,  and  Wimbush  (Ref.  I II -3a, b) . This  chart 
places  the  position  of  the  Northeast  Pacific  M2  amphi- 
drome  at  2 7°N,  135°W.  Schwiderski ' s chart  places  this 
amphidrome  at  26°N,  135°  30'W.  The  figure  shows,  in 
addition  to  the  cotidal  lines,  corange  contours,  which 
are  lines  of  equal  tide  height  with  respect  to  the 
amphidrome.  Both  the  magnitude  and  the  direction  of 
the  tidal  current  vector  are  shown. 


1 1 1-4.  Tidal  Currents  Along  the  Oahu  to  San  Francisco  Path 


Schwiderski  has  calculated  the  M2  tidal  current  components  and  the 
tidal  heights  and  phases  for  the  author  for  four  points  on  the  Oahu  to 
San  Francisco  path  (see  Figure  III -3) . These  are  listed  in  Table  III-4. 
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160  150  140  no  i?o  no 


Figure  III-2.  M2  aotidal  chart,  (from  Munk,  et  al.3  Ref.  III-3a;  repro- 
duced with  permission)  Lines  of  equal  amplitude  H cm  and 
of  equal  Greenwich  epoch  G°  are  shown  by  the  dashed  and 
solid.  curves.  Ellipses  refer  to  the  computed  tidal  cur- 
rents at  the  ellipse  center.  The  heavy  arrow  toward  G = 0 
is  the  tidal  current  when  the  moon  is  over  the  Greenwich 
meridian.  The  dots  indicated  for  30°,  60°,  ...,  give  the 
end  point  of  the  current  vector  for  other  Greenwich  epochs. 
The  dashed  arrow  is  the  current  vector  at  the  time  of  local 
high  tide.  The  inset  shows  the  observed  currents  at  JOSIE 
II.  Tick  marks  on  the  major  axis  of  the  ellipses  corre- 
spond to  current  speeds  at  intervals  of  1 cm/sec. 


III-6 


•Pacific  St.  #1 
37"23.2'N  124®W 

0.64 

303.8 

3.54 

53.9 

52.1 

201.3 

Pacific  St.  #2 
31»56.1'N  139°W 

0.39 

272.7 

2.12 

71.2 

21.5 

257.6 

Pacific  St.  #3 
28°28.2'N  146®W 

0.59 

270.4 

1.81 

72.0 

IS. 9 

298.7 

J 

Pacific  St.  #4 
25»19.2'N  1S2°W 

0. 75 

275.6 

1.71 

7S.4 

16.6 

333.4 

where  (Au>  4u)  = 

Amplitude  and 

phase  of  east 

velocity  u 

i 

(Av>  «v)  - 

Amplitude  and  phase  of  north  velocity  v 

(A^,  6^)  = Amplitude  and  phase  of  tidal  height  £ 

All  phases  4 relative  to  Greenwich  meridian  for  which 
4 * 1*  * 2.0701  min  (time) 

•Marked  velocity  constants  unreliable. 
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IV . BOTTOM  TOPOGRAPHY  OF  THE  ACOUSTIC  PATH 
IV-1.  Effects  of  Bottom  Topography  on  Acoustics 

Topographic  features  may  contribute  to  multipath,  reverberation  and 
reduced  signal  levels.  Topographic  characteristics  relate  to  the  many 
charted  and  uncharted  seamounts  in  the  Pacific.  The  longer  the  acoustic 
path,  the  greater  the  probability  of  multipaths  due  to  reflections  from 
seamounts.  A presentation  of  this  situation  in  the  extreme  is  shown  in 
Figure  IV-1,  taken  from  Sheehy  and  Halley  (Ref.  IV-1).  In  1955  they 
monitored  the  arrivals  from  a nuclear  explosion  in  the  Southeast  Pacific 
several  hundred  miles  southwest  of  San  Diego.  Signals  were  received  at 
three  locations:  Kaneohe,  Pt.  Sur,  and  Pt.  Arena.  The  reflected  signals 
from  islands  and  seamounts  were  recorded  for  2 hours  at  Kaneohe,  and  for 
4 hours  at  Pt.  Sur.  There  were  more  unknown  reflectors  than  there  were 
identifiable  ones. 

IV-2.  Sound  Channel  Signal  Blockage 

Another  effect  of  seamounts  lying  in  the  deep  sound  channel  of  the 
acoustic  path  is  sound  channel  signal  blockage.  Dr.  Henry  Fleming  of 
NRL  furnished  Figure  IV-2  as  an  example,  based  on  actual  measurements 
for  another  path.  The  acoustic  ray  systems  shown  have  been  calculated 
for  two  cases  of  seamount  configuration  in  a 1700  km  path  (not  ours)  for 
a frequency  of  110  Hz.  In  the  case  of  the  upper  figure,  no  signal  was 
received;  in  the  lower  figure  a signal  was  received.  The  positioning  of 
the  transmitter  in  depth  and  range  would  appear  to  be  quite  important 
with  respect  to  the  density  of  rays  clearing  the  seamounts. 

IV-3.  Bottom  Topography  of  the  Acoustic  Path 

The  NORDA  SYNBAPS  program  (Ref.  IV-2)  has  provided  40  great  circle 
radial  bathymetric  profiles  spaced  1°  apart  in  bearing  from  Oahu,  Hawaii, 
to  the  West  Coast  and  approximately  centered  on  the  San  Francisco  radial. 
The  bearings  of  these  radials  are  shown  in  Figure  IV-3. 

SYNBAPS  is  the  acronym  for  "synthetic  bathymetric  profiling  system." 
The  bathymetric  information  is  provided  in  both  machine-plot  and  listing 
form.  The  listing  gives  the  depth  in  meters  at  1-mile  intervals.  A 
plot  is  shown  in  Figure  IV-4  for  the  Oahu  to  San  Francisco  path.  Block- 
age effects  are  not  expected  to  be  severe  for  this  path.  However,  the 
path  skirts  seamount  chains,  and  a slight  change  in  path  orientation 
could  produce  severe  blockage.  If  profiles  are  needed  for  specific 
paths,  these  can  be  provided  within  48  hours. 

Data  can  be  provided  on  magnetic  tape  also.  In  preparing  the  his- 
torical oceanographic  data  of  temperature,  salinity  and  sound  speed, 
bathymetry  for  the  Makapuu  Point  to  San  Francisco  path  will  also  be 
provided  from  the  Scripps  Institution  of  Oceanography  data  file. 

Ref.  IV-3  provides  contours  for  the  North  Pacific  Ocean  and  can  be 
used  also. 


IV-1 


1 


<*b  r«  ARBITRARY 


Figure  IV-2.  Effect  of  seamounts  on  acoustic  signal  transmission:  range  of  1700  km, 
110  Hz.  ( provided  by  H.S.  Fleming,  Naval  Research  Laboratory;  not  the 

path  of  this  study) 


>f  40  SYNBAPS  profile  runs;  locations  are  approximate 
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Figure  IV-4.  SYNBAPS  bathymetric  'profile — Makapuu  Point , Hawaii,  to 
San  Francisco. 


IV-4.  References 

IV-1.  "Measurements  of  the  attenuation  of  low  frequency  underwater 
sound,"  M.J.  Sheehy  and  R.  Halley,  J.  Acoust.  Soc.  Am.  29, 

464-469  (1937) . 

IV-2.  "Synthetic  bathymetric  profiling  system,"  NAVOCFANO  TR233,  February 
1973,  DMA  Supply  Depot,  Philadelphia,  Pennsylvania. 

IV-3.  "Bathymetric  atlas  of  the  North  Pacific  Ocean,"  N.O.  Pub.  1301-2-3, 
1973,  NAVOCFANO,  Washington,  D.C. 

General:  Mr.  Roger  Van  Wyckhouse  of  NORDA  generously  provided  SYNBAPS 
listings  and  plots  for  the  required  radials. 
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V.  SOUND  SPEED  PROVINCES  AND  PC HAN  WATER  MASSES 

V- 1 . Acoustic  Properties  of  the  Oahu  to  San  Francisco  Path 

Figure  V-l  shows,  for  January,  the  average  depth  of  the  mixed  layer  or 
surface  sound  channel,  the  critical  depth,  and  the  depth  of  the  sound  chan- 
nel axis  along  the  path.  These  were  determined  from  Dr.  Eugene  Pudeszwa's 
"Sound  Speed  Profiles  for  the  North  Pacific  Ocean,"  Ref.  V-l.  In  preparing 
this  atlas,  Podeszwa  tried  to  find  the  fewest  number  of  independent  areal 
subdivisions  that  would  furnish  this  information  on  a monthly  basis  with  an 
acceptable  error  in  sound  speed.  For  our  path  the  number  of  independent 
provinces  turns  out  to  be  eight,  corresponding  to  the  eight  sections 
between  the  points  shown  in  the  figure.  The  grouping  of  these  eight 
provinces  among  the  three  water  masses  of  the  path  is  also  shown  in  this 
figure.  Dividing  the  acoustic  path  length,  3835  km,  by  eight,  gives  an 
average  length  of  about  480  km  for  a province.  This  suggests  that  we 
may  be  talking  here  about  stable  mesoscale  features.  Mesoscale  features 
are  discussed  in  Section  VI I — 1 . 

V-2 . Sound-Speed-Profile  Provinces  for  the  Central  and  North  Pacific  Ocean 

Figure  V-2  shows  how  the  acoustic  path  traverses  eight  sound-speed- 
profile  provinces  compiled  for  the  depth  interval  0 to  4500  ft  (Ref.  V-l). 
It  can  also  be  seen  that  the  shapes  of  the  provinces  follow  those  of  the 
current  systems  and  water  masses:  westerlies,  California  Current,  east- 
erly trade  winds,  and  a diagonal  wind  shear  and  current  line  passing 
through  the  Hawaiian  Islands.  Table  V-l  shows  how  the  Podeszwa  sound- 
speed-profile  province  classification  fits  into  the  North  Pacific  water 
mass  classification  using  the  T-S  diagram  classification  of  Figure  1 1 - 2 . 

Table  V-2  lists  three  acoustic  parameters  along  the  path:  mixed 
layer  depth,  deep  sound  channel  axis  depth,  and  critical  depth.  The 
critical  depth  is  taken  as  that  depth  where  the  sound  speed  equals  the 
surface  sound  speed.  It  is  important  in  specifying  convergence  zone 
propagation  for  a near-surface  source.  Note  that  Podeszwa's  tables  do 
not  allow  us  to  describe  the  seasonal  dependence  of  the  sound  channel 
axis  depth.  However,  this  information  may  be  gotten  from  the  summaries 
of  historical  data  obtained  from  the  National  Ocean  Data  Center. 

V - 3 . Summary  of  Historical  Monthly  Sound  Speed  Data 

Lists  of  historical  monthly  sound  speed  data  were  provided  at  stand- 
ard depths  by  the  National  Ocean  Data  Center  (NODC;  NOAA) . The  lists 
were  provided  in  the  2°  x 2°  Modified  Canadian  Square  (MCS)  breakdown  as 
shown  in  Figure  V-3.  The  lists  are  shown  as  Tables  V-3  through  V-21. 

There  is  no  list  for  MCS  1215,  2°  x 2°  #26  because  no  data  existed  in 
the  NODC  files  for  this  region. 
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PODESZWA'S  PROVINCES 

B-18  B- 16  B-14  C-ll  C-IO  D-8  D-7 


-i 1 


JANUARY 

■ PODESZWA  MIXED  LAYER  DEPTH  (m) 

• PODESZWA  SOUND  CHANNEL  AXIS  DEPTH  (m) 
» CRITICAL  DEPTHS 

BOTTOM  TOPOGRAPHY 

HORIZONTAL  TO  VERTICAL  SCALE:  1000:1 


4000h 


4400 1- 1 


»i  l i i ill  i i i l I I I 1 1 1 1 

8 12  16  20  24  28  32  36  40 


DISTANCE  ( I02  km) 

Figure  V-l.  Acoustics  properties  of  path  for  January  from  Makapuu  Point 
Hawaii , to  San  Francisco. 
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Table  V-l.  Relation  between  North  Pacific:  water  mass  class  and  sound- 

speed-profile  provinces  (acoustic  path  between  Makapuu  Point, 
Hawaii,  and  San  Francisco) . 


Distance  Along  Path 


Sound-Speed-Profile 
Province  (Podeszwa) 


Water  Mass  Class 


0 

- 949 

km 

F-4  l, 

Eastern  North  Pacific 

949 

- 1755 

km 

B-18  f 

Central  Water 

1755 

- 2480 

km 

B-16  ) 

Northeast  Pacific 

2480 

- 2889 

km 

B-14 

Transition  Zone 

2889 

- 3210 

km 

C-ll  ) 

3210 

- 3380 

km 

C-10  ) 

California  Current 

3380 

- 3570 

km 

D-8  } 

(modified  subarctic 

3570 

- 3835 

km 

D-7  I 

water) 

Table  V-2.  Seasonal  dependence  of  mixed  layer  depths,  m,  and 

critical  depths,  m,  for  Oahu  to  San  Francisco  path.  * 

Podeszwa ' s Sound  Speed  Provinces 


F-4 

B-18 

B-16 

B-  14 

C-ll 

C-10 

D-8 

D-7 

Distance  Along 

Path  to  Outer 
Boundary,  km 

949  km 

1755  km 

2480  km 

2889  km 

3210  km 

3380  km 

3570  km 

3835  km 

Mixed  Layer 

Depth,  m 

30.5  m 

b 1 . 0 m 

122  m 

107  m 

6 1 . 0 m 

68 . 6 m 

30.5  m 

45.7  m 

Deep  Sound 

Channel  Axis,  m 

807.7 

731.5 

701.0 

t»85.8 

670.6 

685 . 8 

"*31.5 

609 . 6 

Critical  Depth,  m 

4334 

4127 

3596 

3398 

2969 

2895 

274  3 

2194 

Mixed  Layer 

Depth,  m 

7b.  2 

45.7 

107 

61.0 

61.0 

61.0 

22.9 

45.7 

Deep  Sound 

Channel  Axis,  m 

807 . 7 

731.5 

701.0 

685.8 

670.6 

685.8 

731.5 

b09 . 6 

Critical  Depth,  m 

4374 

4267 

3383 

3109 

2774 

2652 

2743 

2179 

Mixed  Layer 

Depth,  m 

15.2 

15.2 

21.3 

15.2 

33.5 

15.2 

18.3 

15.2 

Deep  Sound 

Channel  Axis,  m 

807.7 

731.5 

701.0 

685.8 

670.6 

685.8 

731.5 

609 . 6 

Critical  Depth,  m 

4542 

4587 

3962 

.3901 

3597 

3383 

3185 

2850 

Mixed  Layer 

Depth,  m 

30.5 

33.5 

45.7 

30.5 

38.  1 

24.4 

30.5 

27.4 

Deep  Sound 

Channel  Axis,  m 

807.7 

731.5 

701.0 

685.8 

670.6 

685.8 

731.5 

609 . 6 

Critical  Depth,  m 

4572 

4557 

1 

4191 
| 

4008 

3566 

3459 

3414 

3109 

.January 


/ April 


October 


•Podeszwa 's  tables  do  not  provide  seasonal  dependence  for  the  deep  sound  channel  axis. 
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Figure  V-3.  Modified  Canadian  Square  designators ; 

Oahu  to  San  Francisco , great  circle  path. 


Figure  V-4  shows  a plot  of  the  monthly  values  of  mixed  layer  depths 
and  the  deep  sound  channel  axis  depths,  based  on  the  listings  for  these 
2°  x 2°  squares.  In  this  case  there  were  not  enough  deep  data  to  obtain 
monthly  values  of  critical  depths.  In  plotting  Figure  V-4,  the  distance 
of  the  squares  along  the  path  has  been  taken  nominally  as  the  distance 
of  the  square's  center  from  Makapuu  Point  at  21°19'N,  157°39'W.  These 
distances  are  shown  on  the  tables  of  data.  The  geographical  position 
of  San  Francisco  has  been  taken  as  37°49'N,  122°25'W.  The  great  circle 
length  of  the  entire  path  is  3835  km,  or  2071  n.  mi. 


Dr.  Tom  Davis  of  NAVOCEANO  has  compiled  digital  tape  listings  at 
standard  depths  of  historical  temperature,  salinity  and  sound  speed  data 
for  four  seasons  on  a 30  minute  grid  for  any  point  in  the  oceans  of  the 
Northern  Hemisphere.  He  also  provides  an  error  statement  of  these 
values.  The  beauty  of  the  magnetic  tape  listings  of  sound  speed  is  that 
they  may  be  combined  with  synoptic  XBT  temperature  data  and  surface 
extrapolated  salinity  data  to  provide  an  on-the-spot  sound  speed  profile 
to  the  bottom.  This  procedure's  application  to  this  program  is  discussed 
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DISTANCE  (km) 


0 4 8 12  16  20  24  28  32  36  40 
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Figure  V-4.  Seasonal  changes  in  mixed  layer  depth,  m,  and  deep  sound 
channel  axis  depth,  m,  for  Oahu  to  Scm  Francisco  path; 
NODC  historical  data. 
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Distance  from  Makapuu  Pt . to  center  of  square 

The  maximum  and  minimum  sound  velocities  are  underlined,  showing  the  depths  of  the  mixed 
layer  and  the  sound  channel  axis,  respectively. 


Table  V-4. 

10  Sq  1117  2 Sq  26  (199  km)  Historical  Data,  NODC  Sound  Speed  Median 

(Modified  Canadian)  Profiles,  m/s;  Makapuu  Pt.  to  San  Francisco 
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Table  V-18. 

10  Sq  1214  2 Sq  48  (3175  km)  Historical  Data,  NODC  Sound  Speed  Median 

(Modified  Canadian)  Profiles,  m/s;  Makapuu  Pt.  to  San  Francisco 
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VI.  SAMPLING  OF  TLMPhRATURL  FIELD 

VI-1.  An  Example  of  a Vertical  Temperature  Cross  Section  from  Oahu 
to  San  Francisco  (Sampling  Interval  Along  Path  is  120  km) 

There  are  many  interesting  features  of  this  unique  region  for  which 
oceanic  time  and  length  scales  must  be  considered  in  determining  the 
resolution  of  environmental  detail  that  is  required  to  understand  and 
relate  the  environment  to  acoustic  behavior.  While  it  is  a very  complex 
region,  it  has  been  well  studied.  Ocean  Weather  Ship  (OWS)  NOVLMBLR 
occupies  a permanent  station  halfway  along  the  path  at  30°N  140°W. 

Also,  the  path  lies  along  a major  shipping  lane,  and  since  1966  the 
National  Marine  Fisheries  Service  (NOAA)  has  conducted  an  extensive  XBT 
measurements  program  with  commercial  ships  of  opportunity  (Ref.  VI-1). 

For  about  one  week  of  every  month,  a ship  plies  between  the  two  points, 
making  XBT  temperature  drops  down  to  about  460  m,  along  with  surface 
temperature  and  salinity  measurements  every  4 hours. 

Figure  VI -1  shows  an  example  of  a vertical  cross  section  of  iso- 
therms over  this  path,  taken  from  Fishing  Information  (Ref.  VI -2).  It 
should  be  noted  that  the  closeness  of  an  acoustic  path  to  a shipping 
lane  is  a mixed  blessing  since  there  will  be  more  shipping  noise  to 
contend  with  than  usual.  The  isotherms  shown  are  used  to  study  features 
that  last  much  longer  than  the  5 or  6 days  that  it  took  to  obtain  the 
profiles.  In  the  case  at  hand,  the  sampling  interval  is  about  once 
every  4 hours,  or  every  120  km  if  the  ship  is  moving  at  30  km/h  (16  kn) . 

V I - 2 . A Vertical  Temperature  Cross  Section  from  Oahu  to  San  Francisco 
(Sampling  Interval  of  30  km) 

A more  detailed  plot  (Figure  VI-2)  can  be  made  if  an  XBT  is  dropped 
once  an  hour  for  5-1/2  days.  At  30  km/h  we  get  a sample  spacing  of  30  km. 
Obtaining  XBT  samples  at  this  rate  is  quite  taxing  and  becomes  expensive. 
Also,  at  a ship's  speed  of  50  km/h,  one  is  restricted  to  using  a relatively 
shallow  XBT,  e.g.,  the  T-4  which  goes  down  to  460  m.  It  is  possible  to 
obtain  XBT  data  down  to  760  m (T-7)  by  slowing  down  somewhat  to  27.8  km/h 
(15  kn) . However,  to  obtain  data  down  to  1830  m requires  the  use  of  the 
T-5  XBT  at  a speed  of  11.1  km/li  (6  kn) . The  T-5  is  almost  twice  (1.8) 
as  expensive  as  the  T-7  and  2.2  times  the  cost  of  the  T-4.  Table  VI-1 
lists  the  available  XBT  probes,  and  their  depth  capabilities  and  cost. 

V I - 3 . Accuracy  and  Resolution  of  XBT  Probes 

The  temperature  accuracy  of  the  XBT  probes  is  ±0.2°C  with  a depth 
accuracy  of  ±2%  or  5 m,  whichever  is  greater  (Ref.  VI-3) . The  temperature 
resolution  is  0.1°C,  and  the  time  constant  is  0.5  s.  Thus,  with  a fall 
rate  of  1 m/s  the  probes  respond  to  a change  in  temperature  in  a layer  of 
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Figure  VI- 1.  Honolulu  to  San  Franoisoo , April  1-7,  1976 , HAWAIIAN  QUEEN, 
Voyage  178.  (from  Ref.  VI-2) 
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Figure  VI-2.  Honolulu  to  San  Francisco,  April  27-May  1,  1974,  SS  HAWAIIAN 
ENTERPRISE,  Voyage  89.  (from  Ref.  VI-7) 
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Table  VI- 1.  XBT  probes  (The  Sippioan  Corporation) . 


Model  No. 

Description 

Part  No. 

Unit  Price 
(8/1/78) 

T-4 

Depth  - 460  m (1500  ft) 

Ship  Speed  - 30  kn  (55.59  km/h)* 

207592-1 

$ 25.65 

T-S 

Depth  - 1830  m (6000  ft) 

Ship  Speed  - 6 kn  (11.12  km/h) 

211105-1 

56.85 

T-6 

Depth  - 460  m (1500  ft) 

Ship  Speed  - 15  kn  (24.09  km/h) 

211965-1 

21.30 

T-7 

Depth  - 760  m (2500  ft) 

Ship  Speed  - 15  kn  (27.80  km/h) 

210883-1 

30.90 

T-10 

Depth  - 200  m (660  ft) 

Ship  Speed  - 10  kn  (18.53  km/h) 

213412-1 

15.85 

T-ll 

Finestructure  (FSXBT) 

Depth  - 460  m (1500  ft) 

Ship  Speed  - 6 kn  (11.12  km/h) 

213713-1 

26.10 

XSV-1 

Expendable  Sound 

Velocimeter 

Depth  - 850  m (2790  ft) 

Ship  Speed  - 15  kn  (27.80  km/h) 

213796-1 

66.00 

*1  knot  = 1.853  km/h  = 0.5144  m/s 


water  of  the  order  of  50  cm.  Temperature  error  is  due  mostly  to  servo- 
motor response  in  the  recorder,  the  accuracy  of  the  chart  paper,  and  to 
a smaller  degree,  uncertainties  in  the  thermistor  calibrations  and  the 
wire  resistance  balance.  The  depth  uncertainty  is  due  to  the  individual 
launcher  height,  water  entry  conditions  (surface  wave  height  and  probe 
entry  angle)  and  probe  fall  rate.  These  uncertainties  are  well  under- 
stood, and  attempts  have  been  made  to  limit  them  by  special  calibration 
procedures,  digital  recording,  and  the  like. 

VI-4.  Improvements  to  Standard  XBT  Probe  System 

(a)  XBT  Digital  Recorder  and  Display  System 

Scripps  (SIO)  has  started  to  use  a digital  data  logger  de- 
signed by  Dr.  Mesecar  at  Oregon  State  University  (Ref.  VI-4). 
An  analog-to-digital  converter  interfaces  a PET  microcomputer 
to  a Sippican  launcher  and  probe. 
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C b)  Shipboard  Environmental  Acquisition  System  (SEAS) 

This  system  is  being  developed  and  implemented  by  the  National 
Weather  Service.  It  has  an  XBT  system  with  automated  launching 
devised  by  Plessey  Environmental  Systems.  The  launching  capa- 
bility would  be  1 XBT  probe  every  15  minutes,  with  the  data 
being  recorded  automatically.  Although  the  cost  would  be 
quite  high,  the  sampling  interval  could  be  narrowed,  in  the 
extreme,  to  7.5  km  for  a ship  moving  at  30  km/h. 

(c)  Expendable  Sound  Vclocimeter  (Sippican  XSV) 

A sing-around  sound  velocity  sensor  is  used  to  measure  sound 
speed  to  an  accuracy  of  ±0.25  m/s  down  to  850  m.  Together 
with  a temperature  probe,  this  device  would  be  capable  of 
measuring  density  vs  depth  directly. 

VI -5.  Merging  of  XBT  Data  with  Historical  Deep  Data  Records 

There  are  many  different  procedures  (e.g.,  Refs.  VI -5, b)  for  mer- 
ging depth-limited  XBT  data  with  historical  dee])  ocean  data  records  for 
a particular  location.  Where  less  than  mesoscale  sonar  ranges  (~100  km) 
are  involved,  the  procedural  details  appear  to  have  little  effect.  For 
mesoscale  ranges,  however,  synoptic  variability  of  the  near  surface  iso- 
therm depths  may  be  forced  on  the  historical  data  by  improper  extrapola- 
tion of  shallow  XBTs. 

The  NMFS  program  will  use  a technique  that  merges  XBT  and  extrapo- 
lated surface  salinity  data  together  with  historical,  digital,  deep 
ocean  data  files  supplied  by  NAVQCEANO  to  enable  computation  of  dynamic 
height  (Ref.  II 1-2)  and  sounu  speed  profiles. 

V I - (> . Tempera  t ure  Time  Series  vs  Depth 

It  is  highly  desirable  to  obtain  long-term  time-series  temperature 
variability  data  down  to  depths  below  the  sound  channel  axis  using 
thermistor  arrays  at  a suitable  location  along  the  path.  The  use  of  a 
thermistor  array  (Refs.  VI-8,  9,  10)  could  enable  the  study  of  the  re- 
lationship of  mesoscale  time  changes  in  the  ocean  to  the  propagation 
loss  variability  over  the  whole  path.  It  will  also  give  insight  into 
interpreting  the  depth  dependence  of  temperature  fluctuations  measured 
by  the  shallower  XBT  probes  along  the  path. 

As  an  example  of  the  technique,  deep  ocean  moored  thermistor  string 
buoys  were  used  (Ref.  VI -8)  to  study  the  temperature  and  sound-speed 
time-series  vs  depth  for  the  Eleuthera- Bermuda  acoustic  path  (1250  km). 

In  addition  to  energy  spectra  for  these  fluctuations,  vertical  coherence 
was  studied  for  isotherm  displacements.  Fluctuations  of  various  scales 
were  observed. 
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( a ) Small  scale  fluctuations  of  generally  less  than  1 hour  with 
0.1°C  to  0.2°C  in  amplitude.  These  have  isotherm  displace- 
ments of  10  m at  all  depths.  The  noise  background  was  about 
0.02°C  with  a measured  time  constant  of  2 minutes. 

(b)  Semidiurnal  (tidal)  fluctuations,  with  displacements  of  20  m 
to  50  m could  be  seen  throughout  the  water  column. 

(c)  Eddies  with  3-day  periods,  a 1°  to  2°C  temperature  range,  and 
a displacement  of  50  to  100  m separated  by  intervals  of  12-13 
days. 

(d)  Larger  scale  eddies  were  indicated  but  not  resolvable  by  the 
40-day  record. 
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General:  Much  of  the  information  in  Section  VI  has  been  derived  from 
numerous  personal  communications  with  Dr.  D.R.  McLain,  (NMFS) . 

He  arranged  the  provision  of  the  "Fisheries"  data  and  gave 
the  author  much  assistance  and  consultation  regarding  the 
state  of  knowledge  of  XBT  data  and  techniques.  Dr.  John  P. 

Dugan  (NRL)  provided  pertinent  information  on  the  limitations 
and  use  of  XBT  data.  Mr.  Alvan  Fisher  explained  his  water 
mass  classification  and  merge  technique.  Mr.  K.IV.  Lackie 
(NORDA)  provided  the  information  on  Sippican  XBT  probes  and 
the  LRAPP  experience  with  these  probes. 
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VII.  MESOSCALE  ISOTHERM-DEPTH  VARIABILITY 


VI I - 1 . Power  Spectrum  of  10°C  Isotherm  Depth  Over  Acoustic  Path 

How  closely  spaced  must  XBT  drops  be  to  resolve  the  time  and  length 
scales  of  importance  to  the  acoustic  problem?  We  saw  earlier  in  Sec- 
tion V-l  that,  according  to  Podeszwa's  sound-speed  profile  classifica- 
tion, there  are  eight  independent  provinces  along  the  path.  These  have 
an  average  length  of  about  480  km.  Mesoscale  features  range  from  100  km 
to  upwards  of  1000  km.  If  we  wish  to  model  these  features  for  estimating 
their  acoustic  behavior,  then  the  XBT  spacing  should  be  about  50  km  or 
less.  We  can  then  determine  the  power  spectrum  of  the  variance  of  a 
mesoscale  feature,  e.g.,  10°C  isotherm  depth,  and  find  the  distribution 
of  energy  density  among  the  mesoscale  wavelengths.  In  particular,  the 
time  and  space  behavior  of  the  sound  channel  axis  would  be  important 
for  modeling  purposes. 

Dr.  T.  Davis  has  calculated  for  the  author  the  spatial  power  spec- 
trum (see  Ref.  V 1 1 - 1 ) for  the  10°C  isotherm  depth  for  a densely  sampled 
(42.4  km)  set  of  data  taken  over  the  Oahu  to  San  Francisco  path  by  Drs. 

R.  Bernstein  (Scripps)  and  C.  Collins  (NSF)  in  April  1974  (Ref.  V 1 1 - 2 ) . 
(See  Figure  VI-2.)  The  power  spectrum  (m-/cpkm)  has  been  plotted  against 
wavenumber,  k,  (cpkm)  as  well  as  wavelength  (km)  in  Figure  VII-1.  This 
spectrum  contains  a number  of  bands  for  wavelengths  less  than  1000  km. 
Although  the  spectral  slope  is  flat  from  500  km  to  2000  km,  the  envelope 
falls  off  generally  as  k-2.  The  author  has  examined  analysis  techniques 
for  the  peaks  shown  in  Figure  VII-1.  These  peaks  appear  to  be  real  and 
not  the  result  of  some  artifact  of  the  data  analysis. 

Wilson  and  Dugan  (Ref.  VI I -5)  made  multiship  XBT  drops  in  two  swaths 
south  and  north  of  the  Kuroshio  Extension  in  November  1975  (south)  and 
May  1976  (north).  The  sampling  data  over  the  1800  km  long  southern  swath 
is  shown  in  Figure  VII-2a,b,c.  The  swath  is  216-270  km  wide  and  was  cov- 
ered by  seven  ships.  The  spacing  between  ships  was  about  45  km,  while 
XBTs  were  dropped  about  every  27  km  along  the  track.  Figure  VII -2a  shows 
a plot  of  surface  temperature  for  the  swath  with  its  cold  and  warm  meso- 
scale features--in  fact  there  are  four  cold  intrusions.  Figure  VI I -2b 
shows  the  same  features  in  depth  by  contouring  the  depth  of  the  12°C 
isotherm  which  ranges  between  250  m and  450  m.  Figure  VI I -2c  shows  a 
vertical  temperature  section  for  one  of  the  ships.  Together,  these  fig- 
ures present  a three-dimensional  view  of  mesoscale  features,  showing  the 
relationship  between  surface  and  subsurface  manifestations. 

Figure  VII-3,  taken  from  the  same  study,  shows  the  spectra  for  10°C 
and/or  12°C  isotherm  depths  for  various  latitudes.  The  spectrum  for  the 
less  energetic  southernmost  latitudes,  29-50. 5°N,  resembles  that  of  the 
Oahu-San  Francisco  path  in  both  level  and  slope. 
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Figure  VII-3.  Spectra  of  isotherm  displacement  in  main  thermocline. 

(x)  12°C  isotherm  for  three  ships  29-30. 5°N;  (+)  12°C 
isotherm  for  three  ships  30 . 5-32°N ; (•)  10°C  isotherm 
for  six  ships  36-38°N.  (from  Ref.  VII-3 ; reproduced 
with  permission ) 


Table  VI I- 1 lists  the  power  spectra  computed  for  every  second  month 
for  mid-1976  to  mid-1977  for  the  standard  NMFS  XBT  drop  spacing  of  120  km. 
The  plots  for  these  spectra  are  not  included  here,  but  they  mostly  show 
k_2  to  k*3  dependence  on  wavenumber,  k.  However,  they  do  not  yield  any 
data  for  mesoscale  features  below  240  km.  There  is  also  a larger  mean 
square  error  of  determining  spectral  amplitudes  as  will  be  discussed  in 
the  next  section.  Such  spectra  could  be  useful  if  there  were  more 
associated  information.  However,  the  recommended  solution  is  to  use 
closer  spacing  in  the  vicinity  of  fronts  or  sharp  changes  in  temperature 
and  salinity. 
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Table  VII-1.  Power  spectra  (m  /cpkm)  for  10°C  isotherm  depth, 
Oahu  to  San  Franoisao. 


Wavenumber 

(cpkm) 

Wavelength 

(km) 

Bimonthly  Series  (Year  1976- 

1977)  - 

120  km  Spacing 

8/76 

10/76 

12/76 

2/77 

4/77 

6/77 

Median 

Set 

2.60  x 10-4 

3840 

103 

118 

157 

156 

146 

2.56 

60.68 

3.91 x lO-4 

2560 

94.9 

133 

7.2 

22.7 

19.4 

5.65 

7.42 

5.21  x lO"4 

1920 

28.7 

1.96 

1.2 

10.9 

4.65 

1.51 

4.15 

6. 15  x lO"4 

1 536 

22 . 2 

3.68 

9.87 

6.20 

2.25 

2.39 

3.32 

7.81  x lO-4 

1280 

15.7 

1.37 

6.82 

4 . 36 

1.39 

5.87 

1.71 

9. 11 x lO-4 

1097 

11.8 

1.47 

8.65 

5.19 

2.49 

10.2 

2.65 

1 .04  x 10-3 

960 

7.48 

2.10 

9.50 

2.13 

3.83 

7.16 

2.24 

1 . 17  x 10'3 

853 

4 . 36 

1.53 

4.47 

1.23 

2.67 

2 . 42 

0.772 

1 .30  x 10'3 

768 

3.17 

0.864 

4.24 

3.89 

0.892 

1.99 

0.888 

1 . 43  x 10"3 

698 

2.48 

0.698 

6 . 66 

6.71 

0.722 

4.94 

1 .13 

1.56  x lO-3 
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6.24 

0.607 

6.00 

4.26 

1.81 

10.2 

1.36 
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591 

9.73 
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5.15 

2.07 

2.16 

10.1 

1.91 
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3.90 
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2.40 

5.82 

1.28 

1.95  x lO'-3 
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0.672 

1 . 12 

1.77 

3.58 

5.  17 

4.80 

0.305 
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0.976 

1.87 
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0.810 

2.62 
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2.73  x 10‘ 
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0.0470 

2 . 86  x lO’3 
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0.296 

0.0861 

1.70 
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2.  18 
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0.833 
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1.55 

0.  141 
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0.134 

3.26  x lO-'5 
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0.216 

1.28 
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0.419 
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0.153 

3.39  x 10*-5 
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1.11 

0.200 

0.977 

0 . 334 

0.342 
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0.110 

3.52  x 10‘3 
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0.548 

0.680 

0.601 

0.271 

0.379 

2.85 

0.030 

3.65  x 10'3 
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1.48 

0.426 

0.288 

1.01 

4.40 

0.00943 

3.78  x lO'3 
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0.248 

1.49 

0.204 

0.750 

1.08 

4.43 

0.0131 

3.91  x lO"3 

256 
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0.956 

0.01 84 

1 .17 

0.231 

3 . 59 

0.0594 

4.04  x lO’3 

248 

0.0506 

0.  580 

0.  131 

0.236 

0.103 

2.39 

0.162 

4.17  x lO'3 

240 

0.00465 

0.457 

0.115 

0.0182 

0.390 

4.37 

0.  157 
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VII-2.  Sampling  Error  in  Vertical  Displacement 

Another  necessary  procedure  is  to  find  the  sampling  error  in  ver- 
tical amplitude  (Ref.  VII-1).  This  tests  the  significance  of  the  meas- 
ured isotherm  fluctuation  amplitudes  for  determining  power  spectra. 

From  Table  VII-2  we  can  only  say  that  the  rms  error  in  vertical  ampli- 
tude for  2-data-interval  spacing,  or  85  km,  would  be  about  10  m with 
respect  to  the  error  for  the  42.4  km  spacing.  The  total  rms  error  in 
vertical  amplitude  is  the  root  mean  square  of  the  two  errors.  The 
absolute  error  cannot  be  determined  from  this  data  set. 

Table  VII-2.  rms  sampling  errors  for  10°  isotherm  depths 
(10  degrees  of  freedom) 

XBT  Sample  --  Data  Interval  = 42.4  km 

Sample  Spacing  rms  Error  in  Depth 


42.4  km 
84.8  km 
127.2  km 


0 (assumption) 
9.5  m 
12.2  m 


XBT  Sample  --  Data  Interval  = 120  km 


rms  Error  in  Depth 

0 (assumption) 
16.4  m 

Note:  The  10°C  isotherm  depth  was  read  to  the  near- 
est 2 m in  vertical  amplitude.  This  should 
be  part  of  the  calculated  rms  error.  The  Na- 
tional Marine  Fisheries  Service  should  be 
able  to  supply  values  to  closer  tolerance 
from  the  original  listings. 


For  the  same  set,  Dr.  Saur  (Scripps)  has  found  that  the  amplitude 
of  the  wavelength  oscillations  between  360  and  720  km  is  between  40  and 
80  m in  height.  Those  values  are  significantly  larger  than  the  ~10  m 
error  involved 


Sample  Spacing 

120  km 
240  km 


1 
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VII-3.  Time  Scale  for  Mesoscale  Waves 

If  it  takes  fa  days  to  get  a mesoscale  "snapshot"  of  the  path's 
isotherm  structure,  then  we  may  not  be  able  to  separate  effects  of  time 
scales  of  less  than  12  days.  It  turns  out  that,  theoretically,  baro- 
clinic  Rossby  waves  of  400  km  length  have  a period  of  7.7  months  at  50°N 
latitude.  Bourke  and  Pfeiffer  have  found  nonseasonal  periods  of  7-8 
months  for  cold  water  surges  across  the  Subtropical  Front  from  a study 
of  records  at  OWS  NOVLMBER  (Ref.  VII-4).  The  records  were  taken  between 
1962  and  1970.  Those  authors  believe  that  the  nonseasonal  cold  water 
surges  in  the  upper  250  m at  the  southern  boundary  of  the  Transition 
Zone  are  related  to  nondispersive  baroclinic  Rossby  waves.  For  a dis- 
cussion of  baroclinic  Rossby  waves  in  the  Central  North  Pacific,  see 
Bernstein  and  White  (Ref.  VII-5)  and  Roden  (Ref.  Vll-fa). 

Table  VII-3  summarizes  the  current  estimates  of  time  and  length 
scales  for  mesoscale  waves. 


Table  VII-5.  Typical  time  and  length  scales 
for  mesoscale  waves. 


Horizontal  100  to  1000  km 

Vertical  (10°  isotherm)  -100  m 

Time  7 to  8 months 

(nonseasonal) 


VII-4.  Basin-Scale,  or  Macroscale,  Waves 

It  should  be  noted  that  most  of  the  energy  over  this  path  corre- 
sponds to  an  oceanic  scale  of  the  length  of  the  path  (3835  km)  because 
the  baroclinic  trend  is  monotonic  practically  all  the  way  across  the 
path.  This  can  be  seen  in  Figure  VII-4,  where  median  isotherm  depths 
are  plotted  for  a year  for  10°C,  15°C,  and  20°C.  The  baroclinic  dis- 
placement for  the  10°C  isotherm  is  about  300  m. 


r 

UNIVERSITY  OF  WASHINGTON  • APPLIED  PHYSICS  LABORATORY 


4000 
— I 


Figure  VII-4.  Depth  variability  ranges  of  major  isotherms  for  Makapuu 
Point  to  San  Francisco. 


VII-5.  Required  Spacing  and  Depth  of  XBT  Measurement 

It  would  appear  that  if  we  use  the  NMFS  ships-of-opportunity  pro- 
gram, we  must  try  to  supplement  it  with:  (.1)  XBT  drops  every  hour 
(30  km  spacing) , or  (2)  at  least  use  more  frequent  XBT  drops  in  the 
water  mass  boundary  regions  where  the  largest  changes  occur. 

When  the  riessey  automatic  XBT  launcher  and  recorder  become  avail- 
able, then  we  may  be  able  to  increase  our  dense  coverage  data  in  an 
easier,  but  probably  more  expensive,  manner.  It  would  also  be  helpful  to 
employ  the  1830  m depth  XB'l , since  its  use  does  not  require  much  slower 
ship  speeds  and  since  it  covers  the  important  vertical  regions  of  change 
from  baroclinicity  to  barotropicity  near  the  sound  channel  axis.  Wilson 
and  Dugan  (Ref.  VI I -3)  used  alternate  400  m and  760  m XBTs  in  their 
dense  coverage  of  North  Pacific  swaths.  However,  in  view  of  the  neces- 
sity to  follow  the  displacements  of  the  sound  channel  axis,  the  use  of 
1830  m XBTs  in  some  sequence,  say  one  out  of  three,  along  with  the 
others  is  almost  mandatory.  It  is  important  to  determine  the  inter- 
dependence of  mesoscale  vertical  and  horizontal  displacements,  and  how 
these  affect  the  sound  speed  profile  (Ref.  V1I-7). 

In  Ref.  VII-8,  the  use  of  analytical  modeling  and  limited  data  for  the 
modeling  of  mesoscale  eddy  properties  is  discussed.  In  Chapter  VIII,  we  dis- 
cuss satellite  observations  of  surface  temperature  and/or  dynamic  height  meso- 
scale patterns  as  additional  information  in  relating  the  horizontal  mesoscale 
features  to  the  vertical  displacement  characteristics  in  the  ocean. 
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viii.  satellite  coveral;e  of  an  ocean  acoustic  path 

VIII-1.  Satellite  Oceanography 

The  use  of  satellites  to  monitor  mesoscale  features  of  the  ocean 
surface  on  a continuous  basis  offers  promise  with  regard  to  providing  a 
predictive  capability  for  underwater  sound  transmission  over  very  long 
paths.  The  synoptic  XBT  data  along  paths  like  the  Oahu  to  San  Francisco 
path,  and  an  understanding  of  the  relationship  between  sea  surface  pat- 
terns and  isothermal  or  Lsopycnal  surface  patterns  at  depth,  should 
allow  us  to  extend  our  range  of  coverage  for  modeling  purposes.  The 
required  understanding  can,  of  course,  be  provided  by  historical  data  as 
well  as  theories  and  models  of  ocean  circulation. 

Satellite  areal  coverage  performance  figures  are  given  in  Table 
VIII-1  for  SEASAT-1  (now  i nopcrat i ona 1 ) . 

Table  VIII-1.  SEASAT-1  coverage. 

Pass  Period: 

Spacing  of  Passes: 

Repeat  Coverage: 

Time  to  Cover  Honolulu 
to  San  Francisco  Path 

Relative  Dynamic  Height 
Accuracy : 


90  minutes 
16  km 
5 months 

14-15  days 
±10  cm 


Table  VIII-2,  adapted  from  Reference  VIII-1,  presents  one  estimate 
of  the  required  remote  sensing  accuracy,  resolution,  frequency  of  cover- 
age, and  extent  of  coverage  of  oceanographic  parameters  which  could  lie 
of  importance  in  modeling  extremely  long  range  acoustic  paths. 

Table  VI 1 1 -5,  also  adapted  from  Reference  VI 1 1-1,  provides  an 
abbreviated  list  of  remote  sensing  vehicles,  their  pertinent  instru- 
mentation for  the  parameters  shown,  the  oceanographic  application,  and 
the  sponsoring  agency. 


f 


A useful  data  bank  is  supplied  by  the  GOSST-COMPT  program  of  NOAA- 
NESS  at  Camp  Springs,  MD.  Global  operational  sea  surface  temperatures 
(GOSST)  for  50  km  x 50  km  squares  are  stored  for  data  starting  six  years 
ago.  The  temperature  values  are  five-day  averages. 
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Table  VIII-2.  Needed  parameter's  for  thermohaline  models,  (adapted  from 
Ref.  VIII- 1) 


PARAMETER 

ACCURACY 

RANGE  OF 
VARIABLE 

HORIZ. 
RESOL . 

PERIOD 

AREA  SIZE 

Relative  Sea 

0. 25°C 

10°C 

■HH 

3 h 

200  x 200  km 

Surface  Temperature 

0.5'C 

35°C 

su 

12  h 

Ocean  Basin 

Absolute  Sea 

MEM 

0-35°C 

5 Kin 

3 h 

200  x 200  km 

Surface  Temperature 

m . 

0-35°C 

30  km 

12  h 

Ocean  Basin 

Relative  Surface 

3 cm 

1 m 

■■■ 

200  x 200  km 

Topography 

S cm 

2 m 

££■11 

Ocean  Basin 

Absolute  Wind 

0.2  dynes /cm- 

0-20  dynes/em* 

5 km 

3 h/15° 

200  x 200  km 

Stress  Vector 

0.25  dynes/cm- 

0-20  dynes/cm“ 

30  km 

12  h/45° 

Ocean  Basin 

cm/ sec  (10%) 

4-400 

5 km 

3 h/ 15° 

200  x 200  km 

cm/ sec  (20%) 

4-400 

30  km 

12  h/45° 

Ocean  Basin 

Sea  State 

0-15  in 

5 km 

3 h 

200  x 200  km 

Height 

50  era 

0-15  m 

30  km 

12  h 

Ocean  Basin 

0.05  X. 

0-40  X. 

Salinity 

0.1  X. 

0-40  X. 

30  km 

12  h 

Ocean  Basin 

1-100  mm/h 

mi  ’ 

Rain 

1 mm/h 

1-100  mm/h 

Bill 

Ocean  Basin 

Table  VIII-3.  An  abbreviated  list  of  remote  sensing  vehicles.  (adapted 
from  Ref.  VIII- 1) 


MODE 

DATA 

PARAMETER 

APPLICATION 

SOURCE 

Aircraft 

AXBT 

High  ♦ low  altitude 

B £ W photography 
Laser  profile  tapes 

Rad iomet  ric 
temperatures 

Thermal  profile 
reflectance,  etc. 

Micro-al t itude 

SST 

Acoust i cs 

sea  state  (noise) 

Wave  height  (noise) 
currents,  fronts 

non 

DMSP-5D 

Visual  and  thermal 
photography  and 

CCTs 

Reflectance 
radiom.  temp.  (SST) 

Currents,  fronts, 
eddies 

DN'OM 

NOAA-5 

Visual  and  thermal 
photography  and 

CCTs 

Reflectance 
radiom.  temp.  (SST) 

Currents,  fronts, 
eddies 

NESS 

GEOS- 3 

Alt i meter 

Surface  roughness 

Sat . alt i tude 

Wave  height 
(noi se) 

Sea  surf,  topog. 
(currents,  eddies) 

NASA 

NIMBUS 

Infrared  radiometer 
(11.5  pm  and  0.7  pm) 

SST,  prec.  moist. 

Currents,  fronts, 
etc . 

NASA 
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Such  information  can  be  used  in  conjunction  with  the  monthly  sea 
surface  temperature  data  provided  by  NMFS  for  the  North  Pacific  (Ref. 

VI 1 1 -2).  The  contours  shown  in  Figures  VII I -1,2, 5 are  plotted  from 
data  acquired  and  stored  at  the  Fleet  Numerical  Weather  Central  in 
Monterey,  CA. 

VIII-2.  Mesoscale  Dynamic  Height  Perturbations 

The  spatial  distribution  of  dynamic  height  determines  the  config- 
uration of  the  baroclinic  flow  field.  This  field  consists  of  a mean 
flow  and  perturbations.  The  mean  dynamic  topography  with  respect  to 
1000  decibars  is  well  known  (Ref.  VIII —3)  and  is  shown  in  Figure  V 1 1 1 - 4 . 
The  dominant  features  are  the  subtropical  highs  and  the  subpolar  and 
equatorial  lows,  the  position  and  extent  of  which  vary  with  season  and 
depth.  Typical  length  scales  of  these  features  are  several  thousand 
kilometers.  Table  VI 1 1-4  is  a compilation  of  monthly  data,  provided  by 
NODC  (NOAA)  for  the  Oahu  to  San  Francisco  path.  These  data  are  in  good 
agreement  with  Figure  VI 11-4.  (The  NODC  data  is  referenced  to  1000  m/ 
1500  m;  1 m is  =1  decibar.) 

Roden  (Ref.  VI 1 1-4)  finds  that  data  from  closely  spaced  stations 
about  30  km  apart  along  a 2700  km  path,  with  vertical  sampling  at  3 m 
intervals  down  to  1500  m,  show  perturbations  superposed  on  the  mean 
dynamic  topography.  Wavelike  disturbances  with  length  scales  between 
400  and  600  km  and  heights  of  the  order  of  10  dynamic  centimeters  are 
common  in  latitudes  20°N  to  50°N.  The  wave  amplitudes  are  larger  in  the 
western  than  in  the  central  and  eastern  Pacific.  This  last  observation 
is  in  agreement  with  Section  VII  of  this  report. 

Roden  finds  further  that  the  wavelike  disturbances  extend  down  to 
several  hundred  meters  in  depth  with  an  exponential  decrease  in  wave 
amplitude.  Emery  et  al.  (Ref.  VI 1 1 -5 ) find  significant  mesoscale  per- 
turbations of  the  sound  channel  axis  depth  at  about  800  to  900  m. 

The  baroclinic  currents  associated  with  the  dynamic  height  pertur- 
bations are  an  order  of  magnitude  larger  than  the  mean  currents  and  are 
related  to  the  curl  of  the  wind  stress  field.  One  inference  from  Ref. 

VI 1 1 -3  is  that  a time  series  analysis  of  the  difference  in  sea  level 
readings  corrected  for  tides  between  Oahu  and  San  Francisco  should  yield 
periods  associated  with  the  mesoscale  dynamic  height  perturbations. 

Table  VI 1 1-5  shows  the  contribution  to  dynamic  height  from  1500  m 
to  1000  m in  dyn  cm  for  the  Oahu  to  San  Francisco  path.  There  does  not 
appear  to  be  a trend  by  month  or  distance  for  the  small  amount  of  data. 
The  median  contribution  is  35  dyn  cm. 
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Deviation  of  sea  surface  temperature  (°F)  from  20-year 
mean  (1948-1967 ) , eastexm  North  Pacific  Ocean.  (from 
Ref.  VIII- 2) 
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VIII-3.  Precision  Satellite  Altimetry 

The  author  participated  in  a workshop  on  "Quasi-geostrophic  and 
Mesoscale  Aspects  of  Ocean  Dynamics"  in  connection  with  Ref.  VI 1 1-6, 
during  which  the  required  and  existing  capabilities  of  satellite 
altimetry  were  discussed.  The  results  reported  in  this  section  were 
the  consensus  of  that  meeting. 

Precision  satellite  altimetry  may  allow  the  measurement  of  geo- 
potential topography  over  the  oceans.  The  recent  launching  of  SEASAT-1 
demonstrated  that  a short -pulsed  altimeter  is  capable  of  resolving 
±10  cm  in  change  of  surface  elevation.  This  is  an  order  of  magnitude 
better  than  that  aboard  GEOS-3. 

However,  the  measurement  of  absolute  surface  elevation  requires 
a priori  knowledge  of  the  shape  of  the  geoid  to  within  the  resolutions 
required  for  measuring  surface  elevation.  The  best  geoid  to  date  is 
accurate  to  within  ±1.0  m.  An  accurate  local  geoid  can  be  determined  by 
removing  geoidal  variability  through  data  processing.  The  variability 
of  the  ocean  surface  elevation  with  respect  to  such  a local  geoid  can 
then  be  determined  accurately  to  within  instrumental  errors. 

The  information  in  Table  VI1I-6  is  a condensed  compilation  of  the 
various  ocean  systems  that  have  measurable  geopotential  heights,  and 
indicates  the  required  satellite  altimeter  capability. 

The  working  group  reviewed  the  major  ocean  systems  in  terms  of 
their  geopotential  heights,  horizontal  length  and  time  scales,  and 
current  magnitudes.  The  results  of  the  discussions  are  summarized  in 
the  first  four  columns  of  Table  VI 11-6.  The  required  altimeter  capabil- 
ities for  measuring  the  physical  properties  of  ocean  dynamical  systems  are 
shown  in  the  next  three  columns.  Different  processes  arc  impacted  dif- 
ferently by  geoidal  errors,  tidal  elevation  changes,  orbital  errors,  and 
instrumental  techniques.  These  effects,  investigated  by  other  working 
groups  of  the  Symposium,  are  summarized  in  the  last  four  columns  in  the 
table. 

From  such  considerations  and  the  information  summarized  in  Table 
VI 1 1 - 1 , three  categories  of  oceanic  processes  are  identified:  (a)  those 
that  can  be  studied  with  existing  technology,  (b)  those  that  require  an 
improved  altimeter  resolution  (±5  cm  expected  in  5 years),  and  (c)  those 
that  require  an  improved  best  geoid  to  ±10  cm  (expected  in  10  years). 

(a)  Ocean  Systems  That  Can  Be  Studied  with  Existing  Technology 

1.  Mesoscale  eddies 

2.  Western  boundary  currents 

3.  Antarctic  and  circumpolar  currents 

4.  Curvature  change  in  major  current  systems 

5.  Seamounts  and  trenches 


VI 1 1 -9 
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VIII-1. 

VI 1 1-2 . 

VI 1 1 -3 . 

VI 1 1 -4 . 

VI 1 1-5 . 

VIII-6. 

General : 


UNIVERSITY  OF  WASHINGTON  • APPLIED  PHYSICS  LABORATORY 


(b)  Ocean  Systems  That  Require  ±5  cm  Altimeter  Resolution 

1.  Eastern  boundary  currents 

2.  Upwelling 

3.  Equatorial  currents 

(c)  Ocean  Systems  That  Require  ±5  cm  Altimeter  Resolution  and 
±10  cm  Best  Global  Geoid  (Possible  in  10  Years) 
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IX.  EXTREMELY  LONG  RANGE  ACOUSTIC  PROPAGATION 


I X- 1 . Time-Dependent  Acoustic  Effects 


The  present  fixed  acoustic  path  under  study  between  Oahu  and  San 
Francisco  (3835  km,  133  Hz)  is  the  longest  fixed  path  ever  attempted. 

The  previous  longest  fixed  path  was  from  Eleuthera  to  Bermuda  (1250  km, 
406  Hz).  The  mean  propagation  loss  was  about  120  dB  for  the  1250  km 
path.  For  the  longer  path,  we  estimate  an  additional  loss  of  32  dB 
due  to  cylindrical  spreading,  while  there  is  expected  to  be  8 dB  less 
absorption  loss  for  the  133  Hz  signal  even  though  the  path  is  three 
times  as  long.  Thus  the  expected  propagation  loss  is  144  dB.  The 
travel  time  over  the  longer  path  is  expected  to  be  42.6  minutes.  Pres- 
ently available  acoustic  models  are  expected  to  explain  any  important 
deviations  from  these  values. 

For  the  Eleuthera-Bermuda  path,  the  propagation  loss  and  phase 
record  for  cw  signals  were  recorded  for  a period  of  five  months  (Ref. 

I X- 1 ) . Long-period  phase  and  amplitude  fluctuations  were  observed  to 
be  qualitatively  similar  to  the  short  path  (50  to  75  km)  results  for  the 
Straits  of  Florida.  Phase  fluctuations  associated  with  ocean  tides  were 
a marked  feature  of  the  data.  Shorter  period  fluctuations  of  a few 
minutes  to  a few  hours  were  also  observed. 


Temporal  fluctuation  data  are  relatable  to  environmental  fluctua- 
tion data  and  are  suitable  for  the  testing  and  validation  of  environmental- 
acoustic  propagation.  Models  are  best  obtained  in  fixed-source,  fixed- 
receiver  studies  with  cw  or  periodic  sequence  signals  (Ref.  IX-2).  Rela- 
tionships between  acoustical  and  environmental  fluctuations  have  been 
studied  for  periods  varying  from  a fraction  of  a second  to  a substantial 
fraction  of  a year.  Generally  the  larger  the  time  scale  of  the  environ- 
mental phenomenon,  the  larger  the  amplitude  and  phase  fluctuations  that 
result  from  it. 

Transmission  loss  varies  both  spatially  and  temporally.  Knowledge 
of  spatial  variations  yields  the  maximum  range  of  detection  and  hydro- 
phone array  design.  Knowledge  of  temporal  variations  is  required  to 
improve  signal  processor  design  and  to  predict  system  performance.  Also, 
a study  of  temporal  fluctuations  yields  information  about  the  oceanic 
and  atmospheric  environment. 

Acoustic  amplitude  and  phase  fluctuations  have  been  related  to: 
surface  waves,  internal  waves,  diurnal  and  semidiurnal  tidal  phenomena, 
Rossby  waves,  solar  heating,  changes  in  tidal  phenomena  due  to  lunar 
declination,  movements  of  large  air  masses  and  associated  pressure  dif- 
ferences, seasonal  effects,  currents,  and  f inestructure. 
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Figure  I X- 1 (Ref.  IX-1)  shows  the  two-day  mean  and  variance  of  the 
transmission  loss  over  the  1250  km  El euthera-Bermuda  path  extended  over 
a 5-month  period.  Note  the  long  period  20-dB  change  as  well  as  the 
fairly  stable  7-dB  standard  deviation. 
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Figure  IX-1.  A time -dependent  statistical  characterization  of  the  Bermuda 
amplitude  data  over  a 5-month  period.  Each  data  point  for 
the  mean  and  standard  deviation  and  each  histogram  cover  2 
days  (2880  points)  of  amplitude  data. 


IX- 2.  Range-Dependent  Effects 

There  have  been  numerous  studies  (Ref.  IX-3)  of  the  propagation  of 
sound  through  single  mesoscale  features  such  as  a warm  ring  (eddy) , a 
cold  ring  (eddy),  a front  between  water  masses,  the  Gulf  Stream,  and  a 
Rossby  wave.  Internal  wave  effects  are  not  expected  to  be  important 
(Ref.  IX-4)  , but  this  remains  to  be  seen.  The  mesoscale  features  have 
been  modeled  deterministically  since  changes  occur  in  time  periods  that 
are  long  compared  with  the  acoustic  propagation  time.  The  inverse  prob- 
lem of  describing  the  oceanographic  features  from  the  acoustic  behavior 
is  being  studied  concomitantly  (Ref.  IX-3). 
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When  sound  is  propagated  over  a fixed  path  through  more  than  one 
w iter-mass  type  and  the  associated  fronts,  currents  and  eddies,  there 
are  range-dependent  effects  because  of  changing  bathymetry , tides, 
currents,  and  sound  speed  profiles.  Uniform  attenuation  losses  are 
relatively  small  compared  with  such  effects  as:  convergence  and  diver- 
gence of  rays,  displacements  of  the  sound  channel  axis,  varying  cross- 
sectional  area  of  the  deep  sound  channel,  Doppler  shifts  and  spreading, 
refraction  by  water  currents,  and  blockage  by  bottom  features.  There 
are  also  changes  of  ray  direction  from  sloping  continental  and  island 
shelves.  Redirection  of  surface  scattered  energy  into  the  bottom  at 
high  gracing  angles  is  also  possible  due  to  swells  and  long  gravity 
waves.  Finally,  there  are  certainly  tidal  period  effects.  Proper  mod- 
eling of  long  range  paths  may  require  statistical  treatment  as  well  as 
deterministic  parameters,  as  is  discussed  in  the  next  few  paragraphs. 

Guthrie  et  al.  (Ref.  IX-5)  were  the  first  to  show  the  importance  of 
range-dependent  effects  through  measurements  and  environmental  modeling. 
Measurements  were  made  between  Antigua  and  Newfoundland  at  frequencies 
of  13.89  and  111.1  lie  from  a range  of  100  km  to  a range  of  2800  km. 
Convergence  zones  were  observed  at  long  range,  and  their  positions  were 
found  to  be  frequency  dependent.  The  average  rate  of  increase  of  trans- 
mission loss  with  increasing  range  was  less  than  that  for  cylindrical 
spreading.  Four  oceanographic  regions  were  required  to  model  the  path 
properly.  XBTs  were  taken  every  two  hours,  or  26  km,  to  a depth  of  760  m. 
Of  the  approximately  86  casts,  52  were  used  to  obtain  the  sound  velocity 
profile  as  a function  of  range,  or  one  every  45  km.  The  lengths  of  the 
four  regions  were:  Region  1,  1085  km;  Region  2,  1287  km;  Region  3, 

259  km;  and  Region  4,  10  km.  Bathymetry  was  also  determined  along  the 
track.  It  is  believed  that  for  the  propagation  times  involved,  the 
spectra  plotted  vs  interference  wavenumber  represent  spatial  fluctua- 
tions rather  than  temporal  fluctuations. 

In  Region  2,  1180  to  2467  km,  a near-surface  sound  channel  was 
present  which  did  not  trap  13.89  lit  energy  but  dul  trap  111.1  Hz  energy 
from  the  shallow  source  (21  m) . This  energy  then  traveled  to  long 
ranges  by  various  leakage  mechanisms  into  the  S0FAR  channel.  At  a 
range  of  2790  km,  the  SOFAR  axis  rose  to  100  m,  producing  a marked 
decrease  in  transmission  loss. 

In  another  experiment  from  Midway  Island  (177°W,  28°N)  in  the  North 
Pacific  (Ref.  IX-6),  shots  were  detonated  at  200  m depth  along  a 1650  km 
track  northeast  of  the  hydrophone  located  on  the  SOFAR  axis.  XBTs  were 
cast  every  6 hours,  one-third  of  which  covered  1800  m depth.  Archival 
data  were  used  to  extend  the  1800  m profiles  to  the  ocean  bottom.  Five 
of  these  profiles  spaced  along  the  track  at  0,  609,  962,  1169  and  1651  km 
were  chosen  for  ray  calculations,  because  they  represented  five  distin- 
guishable classes.  ft  was  found  that  multipath  arrival  structure  may  he 
identified  from  ray  calculations  using  the  five  sound  speed  profiles  and 
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bathymetric  data.  The  environmental  changes  affecting  the  arrival  pat- 
tern of  shots  are  of  long  duration.  Signatures  from  shots  more  than  a 
week  apart  combined  to  produce  a highly  regular  pattern.  The  principal 
ray  paths  are  "refracted-refracted"  (RR) . 

In  Ref.  IX-7,  Emery  et  al.  have  modeled  acoustically  one  of  the 
Northwest  Pacific  swaths  of  Wilson  and  Dugan  (Ref.  VII-3).  A range- 
dependent  parabolic  wave  equation  model  was  used  to  calculate  expected 
acoustic  behavior  along  this  path.  At  50  Hz,  the  authors  find  that  the 
resultant  sound  speed  profiles  are  dominated  by  mesoscale  features  which 
are  apparent  as  large  vertical  excursions  of  the  sound  channel  axis  depth. 

A numerical  experiment  was  used  to  model  the  low  frequency  sound  propaga- 
tion along  a 250  km  path.  The  acoustic  intensity  is  plotted  vs  depth 
and  range  for  a hundred  realizations;  that  is,  the  source  is  placed  at 
100  random  points  at  300  m depth.  It  is  found  that  the  standard  devia- 
tion of  the  ensemble  intensity  field  levels  off  to  a saturation  value 
of  5.5  dB  outside  the  main  beam  in  several  tens  of  kilometers.  This 
variability  is  estimated  to  have  length  scales  on  the  order  of  10  km 
in  range  and  several  hundred  meters  in  depth. 

The  Naval  Oceanographic  Office  surveyed  the  same  general  region  a 
year  or  two  later  (July  and  August  1977;  Ref.  IX-8).  They  used  cw 
signals  (89  Hz  and  107  Hz  at  95  m depth)  and  SUS  charges  (90  m and  18  m 
depths)  to  study  the  acoustic  path  which  was  about  1722  km  long.  This 
path  extended  through  the  middle  of  the  swath  region  in  a NW  to  SE 
direction.  The  environmental  objectives  of  this  survey  were:  (1)  to 
determine  the  thermohaline  and  sound  velocity  structures  of  the  Subarctic- 
Subtropical  Transition  Zone,  and  (2)  to  characterize  the  fronts  marking 
this  region. 

Depth  profiles  of  temperature,  salinity  and  sound  velocity  were  ob- 
tained at  26  sites  along  the  path,  spaced  a distance  of  66  km  apart. 

XBT's  were  dropped  hourly  when  strong  horizontal  gradients  were  found. 
Airborne  thermal  surveys  were  also  conducted  to  observe  spatial  and 
temporal  variations  within  the  Transition  Zone.  Bathymetry  along  the 
track  was  obtained  directly. 

Figure  I X- 2 shows  the  existence  of  a warm  core  eddy  about  300  km  in 
diameter  occupying  much  of  this  part  of  the  Transition  Zone  which  has  an 
average  width  of  500-1000  km  across  the  Pacific.  Figure  TX-3  shows  the 
associated  sound  speed  structure. 

For  comparison  with  measurements,  results  are  available  from  range- 
dependent  model  computations  of  the  acoustic  field.  Twelve  profiles  were 
used  to  represent  the  sound  speed  environment  for  this  1600  km  (870  n.  mi.) 
path.  Seven  profiles,  26  km  apart,  were  used  to  cover  the  subarctic  front 
which  was  about  185  km  wide.  Two  profiles  covered  the  subtropical  frontal 
boundary,  while  one  or  two  profiles  covered  the  more  homogeneous  regions. 
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Figure  IX-2.  Vertical  distribution  of  temperature  (°C)  bettieen  35°50'N,  170°04’E 

and  48° 30  'N,  161°53  'E—  22-27  July  1977 ; 1 n.  mi.  = 1.85325  km.  ( from 
Refs.  IX- 8,  11) 
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In  general,  when  a front  is  traversed  from  a warm  to  cold  water  mass,  the 
deep  sound  channel  migrates  upward  along  with  a narrowing  of  the  channel 
cross  section  which  concentrates  the  acoustic  energy.  A front  may  intro- 
duce additional  axially  propagating  paths  from  the  refracted-surface  re- 
flected paths  of  convergence  zones. 

The  Parabolic  Wave  Equation  (PE)  model  (Ref.  IX-9)  and  the  Multiple 
Profile  Program  (MPP)  (Ref.  IX- 10)  have  been  used  to  study  the  expected 
effects  of  the  fronts  and  eddy  of  this  experiment.  Consideration  of  the 
range  dependence  of  the  environment  introduces  a gain  of  about  10  dB 
over  the  constant  environment  for  the  89  Hz  source  at  95  m depth.  (See 
Fig.  IX-4a.)  Although  not  shown  in  this  figure,  the  source  depth  of  95  m 
reduces  the  propagation  loss  over  that  at  18  m by  about  5 dB  in  the  frontal 
region.  An  increase  in  frequency  from  89  Hz  to  107  Hz  increases  the  loss 
by  3 dB. 

Figure  IX-4b  is  a plot  of  the  measured  propagation  loss  data.  It 
may  be  seen  that  at  375  n.  mi.  the  front  produces  a gain  of  about  5 dB 
instead  of  the  calculated  10  dB.  Mote  that  the  transmission  loss  scales 
for  Figures  IX-4a  and  b are  different. 
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Figure  IX-4.  a — calculated  transmission  loss  versus  range — environmental 
effects.  Range  dependent  environment  versus  constant  en- 
vironment. (from  Ref.  IX-8)  b — measured  transmission 
loss  versus  range.  Range  dependent  environment  for 
conditions  specified  in  legend  above.  (supplied  by 
Mr.  R.  Merri field  of  NAV0CEAN0) 
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IX-3.  Attenuation  Coefficients  at  Very  Low  Frequencies 

There  is  inconsistency  in  the  attenuation  data  at  frequencies  below 
100  Hz  (Ref.  IX- 12).  Some  of  the  measured  attenuation  coefficients  exhibit 
a regional  dependence.  One  hypothesis  for  very  low  frequency  attenuation 
is  the  frequency-independent  ray  diffusion  loss  model  due  to  volume  scat- 
tering (Refs.  I X- 13,  14).  Leakage  into  the  bottom  out  of  the  deep  sound 
channel  is  another  hypothesis  (Ref.  I X- 1 5) . This  would  require  a minimum 
at  about  50  Hz  or  less.  The  existence  of  a minimum  attenuation  coeffi- 
cient at  50  Hz  has  also  been  attributed  to  modal  leakage  due  to  near- 
surface scattering  effects  (Ref.  IX-16).  There  is  a strong  suggestion 
that  the  scattering  agent  is  depth  dependent  and  dominant  in  the  upper 
few  hundred  meters  of  the  ocean.  This  scattering  agent  appears  to  be 
related  to  turbulence,  water  velocity  shears,  f inestructure,  and  micro- 
structure. 

Table  IX-1  summarizes  the  magnitude  of  these  regional  attenuation 
coefficients  (adapted  from  Ref.  IX-12). 


Table  IX-1.  Regional  attenuation  coefficients,  dB/km;  f < 100  Hz. 
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The  absorption  coefficient  due  to  the  boric  acid  relaxation  at 
1 kHz  is  regionally  dependent  because  of  pH  variability.  (See  Ref. 
IX-17.)  A typical  pH  and  sound  speed  profile  for  the  Northeast  Pacific 
is  shewn  in  Figure  IX-5.  A relatively  simple  expression  for  estimating 
this  boric  acid  absorption  is  given  in  Ref.  I X- 1 8 along  with  a geo- 
graphic plot  (Fig.  IX-6)  of  associated  contours  for  the  parameter  A 
(f  in  kHz) . 
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Figure  IX-5.  Pacific  Ocean  (Modified  Canadian  Square  1217 ; 30°  to  40°N, 

150°  to  160°W).  A = sound  speed  versus  depth ; B = pH  versus 
depth. 


Figure  IX-6. 

Contour  plot  of  the 
value  A in  the  Pacific 
Ocean  for  use  with  sound 
absorption  equation.  A 
is  determined  from  the 
pH  at  the  SOFAR  axis, 
(from  Ref.  IX- 18 ; repro- 
duced with  permission) 
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IX-4.  Ambient  Noise  Field  at  Very  Low  Frequencies 

Some  data  on  ambient  noise  levels  as  a function  of  depth  for  the 
Northeast  Pacific  Ocean  have  been  published  (Refs.  IX- 19,  20)  from  the 
NORDA-LRAPP  program.  These  data  were  obtained  at  a site  (Site  A, 

30°N,  145°W)  midway  between  Oahu  and  San  Francisco,  and  also  at  a site 
(Site  C,  40°N,  145°W)  about  1000  km  directly  to  the  north.  Considering 
that  the  low-frequency  noise  arose  from  shipping,  the  data  were  used  to 
estimate  the  attenuation  coefficients  for  this  area. 

The  data  were  obtained  during  two  separate  24-hr  periods  and  were 
analyzed  in  1/3-octave  bands  over  the  frequency  interval  10-300  Hz. 

Median  spectral  levels  were  calculated  for  10  sec  intervals.  Figure  IX-7 
shows  the  sound  axis  spectral  noise  levels  for  a 24-hr  period  (period  2) 
at  both  sites  and  at  two  depths,  the  sound  channel  axis  and  the  critical 
depth.  The  results  indicate  a dependence  on  frequency,  hydrophone  depth, 
and  sound  speed  profile. 


Figure  IX-7. 

Median  spectra  for  axial  and 
critical  depths  at  sites  A and 
C for  period  2.  (from  Ref. 

IX- 19 ; reproduced  with  per- 
mission) 
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Figure  IX-8  shows  the  bathymetry  and  sound  speed  structure  along 
the  145°tV  meridian  vs  latitude.  It  can  be  seen  that  the  sound  channel 
axis,  the  critical  depth,  and  the  bathymetry  all  shoal  toward  the  Alaskan 
shelf.  This  figure,  in  conjunction  with  Figure  V-l,  gives  more  perspec- 
tive to  the  problem  of  predicting  acoustic  properties  for  extremely  long 
range  paths  in  the  Northeast  Pacific. 


DEGREES  NORTH  LATITUDE 

30  35  40  45  50  55  60 


Figure  IX-8.  Bathymetric  and  sound-velocity  structure;  14S°W  longitude, 
(from  Ref.  IX- 19 ; reproduced  uith  permission) 


Figure  IX-9  shows  that  the  values  of  the  median  spectral  noise  levels 
at  the  sound  channel  axis  are  consistently  higher  at  the  northern  site. 

The  higher  northerly  noise  levels  are  attributed  to  the  density  of  dis- 
tant shipping.  The  important  point  is  made  that  there  is  a decided 
difference  between  the  axial  noise  levels  and  those  at  critical  depth. 
These  data  are  shown  in  Figure  IX-9  for  both  sites,  and  may  be  used  in 
two  ways  pertinent  to  our  present  objective: 


(a)  The  differences  in  the  two  sets  of  differences 

(4j  axial  - critical)  are  used  to  estimate  the  atten- 
uation coefficient  as  a function  of  frequency  for  the 
area.  These  are  summarized  in  Figure  I X - 1 0 . The  symbol 
aos  stands  for  the  scattering  or  diffusion  attenuation 
coefficient;  aA  is  the  absorption  coefficient;  oqy  is  the 
sum  of  the  two.  The  3 dB  difference  at  130  Hz  for  1000  km 
leads  to  an  estimated  attenuation  loss  of  11  dB  over  a 
3800  km  path. 
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Figure  IX- 9. 

The  difference  between  the 
axial  noise  excess  ( over  the 
critical  depth  value)  at  the 
two  sites  A and  C,  as  a func- 
tion of  frequency.  (from  Ref. 
IX- 19 ; reproduced  with  permis- 
sion) 


i 


I 


Figure  IX-10. 

Low-frequency  attenuation 
coefficients  for  the  North 
Pacific  Central  Water  between 
sites  A and  C.  (from  Ref. 
IX-19 ; reproduced  with  per- 
mission) 
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IX-5 
I X- 1 

I X- 2 

IX-3 


I X-  4 . 


(b)  If  the  3 dB  difference  in  noise  levels  between  the  sound 
axis  and  critical  depth  at  90  Hz  for  Site  A were  to  hold 
everywhere  along  the  path  for  a signal  projected  from 
Oahu,  then  the  shallow  array  at  the  sound  channel  axis 
off  San  Francisco  would  be  at  least  that  much  better  than 
a deep  offshore  array  below  the  critical  depth. 
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X.  CONCLUSIONS 

1.  It  is  feasible  to  use  the  National  Marine  Fisheries  Service  (NMFS, 
NOAA)  ship-of-opportunity  XBT  data  merged  with  the  NAVOCF.ANO  taped 
historical  oceanographic  data  file  to  obtain  sound  speed  profiles 
for  use  in  propagation  calculations  related  to  the  Northeast 
Pacific  fixed-source  fixed-receiver  propagation  study. 

2.  The  shallow  monthly  XBT  runs  must  be  augmented  with  760  m and  1830  m 
runs. 

3.  In  regions  of  large  horizontal  gradients  of  temperature  and/or 
salinity,  XBT  drops  should  be  made  every  30  km  instead  of  every 
120  km.’ 

4.  The  space-time  scale  of  mesoscale  processes  must  be  determined 

in  the  regions  where  propagation  measurements  are  being  performed. 

5.  The  position  and  scale  of  surface  mesoscale  features  may  be  esti- 
mated by  using  NMFS'  surface  temperature  deviation  contours  with 
respect  to  20-year  means. 

6.  Additional  synoptic  XBT  coverage  for  other  locations  in  the  North- 
east Pacific  may  be  obtained  from  the  Fleet  Numerical  heather 
Central's  data  network. 

7.  Long  time-series  of  temperature  at  various  depths  down  to  1500  m 
would  be  useful  in  relating  environmental  variability  to  acoustic 
variabi 1 itv . 

8.  Existing  bathymetric  (SYNBAPS)  and  tidal  (M2)  (Dahlgren)  data  banks 
and  retrieval  systems  appear  to  be  feasible  for  acoustic  modeling. 

9.  Dynamic  height  measurements  obtained  from  satellites  with  +10  cm 
altimetry  capability  may  be  used  to  define  mesoscale  surface  fea- 
tures. 

10.  The  suitability  of  existing  acoustic  models  for  propagation  loss 
and  time-of-arrival  computations  must  be  tested. 
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Attn:  Library 

1 

Crystal  Plaza  *5 

Arlington,  Virginia  20360 

Conmanding  Officer 

Attn:  Dr.  James  Probus,  NAVWAT  03L 

1 

Naval  Intelligence  Support  Center 

4301  Suit  land  Road 

Commander 

Washington,  D.C.  20390 

Naval  Ocean  Systems  Center 

Attn:  Dr.  Johann  Martinek 

1 

Department  of  the  Navy 

Mr.  E.  Bissett 

1 

San  Diego,  California  92132 

Attn:  Dr.  Daniel  Andrews 

1 

STIC  34 

1 

Dr.  Dean  Hanna 

1 

Director 

Mr.  Henry  Aurand 

1 

Applied  Physics  Laboratory 

Dr.  Harry  A.  Schenck 

1 

University  of  Washington 

1013  N.E.  40th  Street 

Commander 

Seattle,  Washington  98105 

Naval  Surface  Weapons  Center 

Attn:  Dr.  T.E.  Ewart 

1 

Acoustics  Division 

Silver  Spring,  Maryland  20910 

Dr.  M.  Schulkin 

1 

Attn:  Dr.  Zaka  Slawsky 

1 

Applied  Research  Laboratories 
University  of  Texas  at  Austin 

Cowander 

P.0.  Box  8029 

Naval  Surface  Weapons  Center 

10000  FM  Road  1325 

Science  and  Mathematics  Research 

Austin,  Texas  78712 

Group  (K05) 

Attn:  Dr.  Lloyd  Hampton 

1 

Dahlgren,  Virginia  22448 

Attn:  Dr.  E.W.  Schwiderski 

1 

Dr.  Charles  Wood 

1 

Dr.  Melvin  J.  Jacobson 

Commanding  Officer 

Rensselaer  Polytechnic  Institute 

Naval  Underwater  Systems  Center 

New  London  Laboratory 

Troy,  New  York  12181 

1 

New  London,  Connecticut  06320 

Dr.  T.G.  Birdsall 

Attn:  Dr.  William  Von  Winkle 

1 

Cooley  Electronics  Laboratory 

Or.  A.  Nuttall 

1 

University  of  Michigan 

Mr.  A.  Ellinthorpe 

1 

Ann  Arbor,  Michigan  48105 

1 

Dr.  D.M.  Viccione 

Mr.  A.  Donn  Cobb 

Nr.  R.  Hasse 

1 

1 
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Dr.  M.A.  Basin 
S . D . P . Inc. 

15250  Ventura  Boulevard 
Suite  518 

Sherman  Oaks,  California  91403 

University  of  Miami 
Rosenstiel  School  of  Marine  5 
Atmospheric  Sciences 
4600  Rickenbacker  Causeway 
Miami,  Florida  33149 
Attn:  Dr.  H.  DeFerrari 
Dr.  W.  Duing 

Dr.  David  Middleton 
127  East  91st  Street 
New  York,  New  York  10028 

Dr.  Donald  Tufts 
Uni  ersity  of  Rhode  Island 
Department  of  Electrical  Engineering 
Wakefield,  Rhode  Island  02881 

Dr.  Loren  Nolte 
Duke  University 

Department  of  Electrical  Engineering 
Durham,  North  Carolina  27706 

Mr.  S.W.  Autrey 
Hughes  Aircraft  Company 
P.0.  Box  3310 

Fullerton,  California  92634 

Dr.  Thomas  W.  Ellis 
Texas  Instruments,  Inc. 

13500  North  Central  Expressway 
Dallas,  Texas  75231 

Institute  for  Acoustical  Research 
Miami  Division  of  the  Palisades 
Geophysical  Institute 
615  S.W.  2nd  Avenue 
Miami,  Florida  33130 
Attn:  Mr.  M.  Kronengold 
Dr.  J.  Clark 
Dr.  S.  Adams 

Mr.  Carl  Hartdegen 
Palisades  Sofar  Station 
Bermuda  Division  of  Palisades 
Geophysical  Institute 
FPO  New  York  09560 

Mr.  Beaumont  Buck 

Polar  Research  Laboratory 

123  Santa  Barbara  Avenue 

Santa  Barbara,  California  93101 

Dr.  M.  Weinstein 
Underwater  Systems,  Inc. 

8121  Georgia  Avenue 

Silver  Spring,  Maryland  20910 

Dr.  Thomas  G.  Kincaid 
General  Electric  Company 
P.0.  Box  1088 

Schenectady,  New  York  12301 

Dr.  C.N.K.  Mooers 
University  of  Delaware 
Newark,  Delaware  19711 


Woods  Hole  Oceanographic  Institution 
Woods  Hole,  Massachusetts  02543 
Attn:  Dr.  Paul  McElroy 
Dr.  R.  Spindel 

Dr.  John  Bouyoucos 
Hydroacoustics,  Inc. 

321  Northland  Averue 
P.0.  Box  3818 

Rochester,  New  York  14610 

Dr.  Richard  James 
c/o  Fleet  Weather  Facility 
4301  Suitland  Road 
Washington,  D.C.  20390 

Atlantic  Oceanographic  5 Meteorological 
Laboratories 

15  Rickenbacker  Causeway 
Miami,  Florida  33149 
Attn:  Dr.  John  Proni 

Research  Triangle  Institute 
Research  Triangle  Park 
Durham,  North  Carolina  27709 
Attn:  Dr.  E.G.  Baxa,  Jr. 

Science  Applications,  Inc. 

8400  Westpark  Drive 
McLean,  Virginia  22102 
Attn:  Ms.  Angela  D'Amico 
Ms.  Lcrna  B lumen 

Bell  Telephone  Laboratories 
1 Whippany  Road 
Whippany,  New  Jersey  07981 
Attn:  Dr.  Bruce  Bogart 
Dr.  Peter  Hirsch 

General  Electric  Company 
Heavy  Military  Electronic  Systems 
Syracuse,  New  York  13201 
Attn:  Mr.  Don  Winfield 
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